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Amphibians occupy a range of habitats throughout the world, but 
species richness is greatest in regions with moist, warm climates. We 
modeled the statistical relations of anuran and urodele species richness with 
mean annual climate for the conterminous United States, and compared the 
strength of these relations at national and regional levels. Model variables 
were calculated for county and subcounty mapping units, and included 
40-year (1960-1999) annual mean and mean annual clim; statistics, 
mapping unit average elevation, mapping unit land area, and estimates of 
anuran and urodele species richness. Climate data were derived from more 
than 7,500 first-order and cooperative meteorological stations and were 
interpolated to the mapping units using multiple linear regression models. 
Anuran and urodele species richness were calculated from the United States 
Geological Survey's Amphibian Research and Monitoring Initiative (ARMI) 
National Atlas for Amphibian Distributions. The national multivariate linear 
regression (MLR) model of anuran species richness had an adjusted coeffi- 
cient of determination (R2) value of 0.64 and the national MLR model for 
urodele species richness had an R? value of 0.45. Stratifying the United 
States by coarse-resolution ecological regions provided models for anurans 
that ranged in R? values from 0.15 to 0.78. Regional models for urodeles 
had R? values ranging from 0.27 to 0.74. in general, regional models for 
anurans were more strongly influenced by temperature variables, whereas 
precipitation variables had a larger influence on urodele models. 


INTRODUCTION 


Amphibian populations appear to be declining worldwide (Houi AHAN et al., 2000, 
Carty ct al, 2001: YounG et al., 2004) A number of possible causes of dechne have been 
proposed including changes in climate {e g . Pot NDS & Cru mp, 1994; Donner ty & CRUMP, 
1998, Por nps et al. 1999), mereased UV radiation (eg. Azrorn & RICHARDS, 1999, 
BLAUSTEIN et al. 2003), habitut loss/fragmentation/alierauon (eg. FalRiG et af. 1995, 
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DEMAYNADIER & Hunter, 1998; Krzysik, 1998; COLLINS & STORFER, 2003), introduction of 
nonindigenous competitive species {€ g., HAYES & JENNINGS, 1986, ROSEN & SCHWALBE, 1995, 
FISHER & SHAFFER, 1996; Kiesecker & BLAUSTEIN, 1998; LAWLER et al., 1999), occurrence of 
contaminants (e.g., BERRILL et al., 1993, BONIN et al., 1997; Davipson et al., 2002, HAYES et 
al., 2002), exposure to pathogens (€ g., LAURANCE et al., 1996, KIESECKER & BLAUSTEIN, 1999) 
and over-harvesting (KOONTZ, 1992: LANNOO, 1996} Many herpetologists believe that com- 
binations of stresses are being placed on amphibian populations (GREEN, 1997; BRITSON & 
THRELKELD, 2000; COLLINS & STORFER, 2003; LANNOO et al . 2003; LITTLE et al , 2003). 


There is widespread acknowledgment that the global climate 1s changing (HOUGHTON et 
al , 2001). Changes in land cover may also affect climate by altering the physical properties of 
the land surface (HAYDEN, 1998; PiFLKE et al , 1999, 2002) Short- and long-term changes in 
chmate have the potential to affect the ranges of individual amphibian species and hence 
species richness in any given locality (THomas et al., 2004) Climatic conditions not only 
directly stress amphibian populations (Dopp, 1997, Pounps et al., 1999, Corn & MUTHS, 
2002), but also may influence their resistance to disease or their ability to withstand attacks by 
environmental pathogens (CAREY & ALEXANDER, 2003). Water availability, air temperature 
and relative humidity can influence amphibian breeding, development, foraging, mobility, 
calling, immune response and habitat availability (DONNELLY & Crump, 1998, GiBBs & 
BREISCH, 2001). Climate also can influence the spread of amphibian pathogens (Daszak et al, 
2003; JOHNSON & CHASE, 2004). 


Amphibians have a long (~350 million years) history of survival under extremes in global 
climate (Carey & ALEXANDER, 2003), yet their life histories (DUFLLMAN, 19994) suggest that 
imdividual amphibian populations may be vulnerable to short-term variations n climate. 
Amphibians occupy a range of habitats throughout the world, but species richness 1s greatest 
in regions with moist, warm chmates (DUELLMAN, 19996) Other natural factors that can 
affect amphibian species richness mclude historical lineages, barriers to migration, interspe- 
cies competition and the avatlability of food, shelter and breeding sites. 


OBJECTIVES AND SCOPE 


This research assesses the degree to which average climatic conditions in the contermi- 
nous United States over the last four decades explain historical patterns of amphibian species 
richness. The primary objectives of the research were to model the statistical relations of 
anuran and urodele species richness with mean annual climate for the contermmous United 
States, and to compare the strength of these relations at national and regional levels, Trends 
in climatic conditions during this period were also evaluated to determine sf they might be 
leading towards more stressful conditions for amphibians te g , decreases in available breeding 
habitat, shortening of breeding season). 

There were limitations or biases implicit in the datasets used for these analyses. First, the 
species occurrence records incorporated into the ARMI National Atlas for Amphibian 
Distributions were not associated with an explicit time period and are best described as an 
histone compilation of ovcurrence records. Second. the species richness estimates were 
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compared to climate statistics averaged for 1960-1999. Climate from this time period may not 
match longer-term averages or averages from other time periods. Third, the species richness 
estimates [http //www.mp2-pwre usgs.gov/armiatlas/] were based on mapping units 
{counties/subcounties} that are not of uniform size, resulting m the potential for an inherent 
bias towards a larger number of species occurring in larger counties. Fourth, neither the 
weather stations nor the spatial variability of weather were uniformly distributed across the 
United States, so the quality of information varied from mapping unit to mapping unit Fifth, 
the mean elevation was used for each mapping unit, but the concept of average elevation may 
not be useful in mountainous areas. No attempt was made to account for the patchiness of the 
distribution of species within a county/subcounty (KuEsTER, 1971), because such data do not 
exist for most of the United States, and no attempt was made to account for the effect of 
climate variation on species that prey upon amphibians or that amphibians consume, Despite 
these hmitations, and because strong relationships between climate and ecosystem develop- 
ment are widely recognized (WALTER, 1973; FORMAN & Gopron, 1986; MONSERLD & 
Leemans, 1992), it was appropriate to expect that relations between climate and amphibian 
species richness would emerge from the analysis. 


METHODS 


SOURCE AND PROCESSING OF AMPHIBIAN DATA 


Species richness estimates for anurans and urodeles were derived from the ARMI 
National Atlas for Amphibian Distributions (hereafter called “atlas”) [http.//www.mp2- 
pwre usgs.gov/armiatlas/], which uses a combination of counties and subcounties as a spatial 
framework for documenting the geographic occurrence of the nearly 300 species of amphi- 
bians currently recognized in the United States (Lannoo et al., 2005). Counties are used as 
mapping units for all but five Western states (Arizona, California, Nevada, Oregon, and 
Washington), for which subcounties are used to help overcome the wide disparity in county 
sizes across the nation. The atlas ıs a compilation of both current and historic records of 
amphibian occurrences, bounded by no explicit time period. The records are from peer- 
reviewed scientific hterature, museum vouchers, state and regional herpetological atlases, and 
other confirmed and validated observations. Data sources vary by state and are not standard- 
azed in their geographic precision Thus, some records in the atlas may represent assumed 
presence, as from a range map. whereas other records represent vouchered specimens with 
specific location information. Because the atlas database incorporates a county/subcounty 
coding system that follows Federal Information Processing Standards, a geographic informa- 
tion system (GIS) was used to link species occurrence records with a digital map of county 
and subcounty polygons [http//www census sov/geo/www/cob/scale html]. Species richness 
was calculated for anurans and urodcles by tallying the number of species recorded as 
occurring within each map unit. 
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SOURCE AND PROCESSING OF CLIMATE DATA 


Estimates of 1960-1999 mean annual and annual mean climate statistics were calculated 
from approximately 7,500 National Weather Service first-order and cooperative temperature 
stations, and 11,500 National Weather Service first-order and cooperative precipitation 
stations. First-order stations are operated by professional staff and report a comprehensive 
array of weather variables each hour. Cooperative sites are more numerous, but generally only 
make once-daily observations of a few weather variables (e g , minimum and maximum dally 
temperature and precipitation) These data were extracted from the National Climate Data 
Center Summary of the Day Dataset and have been quality controlled by the National 
Chmate Data Center (EIsCHEID et al., 2000; CLARK et al , 2004). Estimates of mean annual 
and annual mean climate statistics (tab. 1) for each county/subcounty were calculated using 
multiple linear regression (MLR) models. The MLR method was used to distribute the 
climate statistics (dependent variable) calculated at each station to each county/subcounty 
based on the “XYZ” value (longitude X, latitude Y, and mean elevation Z, respectively) of the 
county/subcounty polygon centroid (Hay et al., 2000, Hay & Mccase, 2002; Hay & CLARK, 
2003). The MLR equation [1] was developed for each dependent variable (climate statistic, 
“CS”) using the independent XYZ variables from a set of National Weather Service climate 
stations, 


CS =b; x +bzy +b; z+ b [1] 


The MLR equations were computed to determine the regression surface that described 
the spatial relations between the dependent CS and the independent XYZ variables, Equation 
[1] describes a plane ın three-dimensional space with slopes b,, b, and b, intersecting the CS 
axis at bg, The best MLR equation for each CS did not always include all the independent 
variables, 


To estimate the chmate statistics for each county/subcounty (CNTY), the following 
procedures were followed: first, mean daily CS and corresponding mean XYZ values from a 
set of stations (STAMEAN) were used with the slopes of the MLR from equation [1] to 
estimate a unique y-intercept (byest, see equation [2]), and second, equation (3] was solved 
using the coefficients (b,, b, and b,) from equation [1], byest from equation [2], and the XYZ 
values of the CNTY. 


byest = CS(STAMEAN) (b; x(STAMEAN) + b, YSTAMEAN) + b, 2(STAMEAN)) [2] 
CS(CNTY) = byest + b, (CNTY) + b, YCNTY) + b, ACNTY) [3] 


The set of stations comprising the STAMEAN in each calculation were chosen from the 
20 closest stations to the CNTY. Outliers (1.e., stations determined to be too far away from the 
data site or residing in another physiographic region) were not used ın the STAMEAN 
calculation The same MLR equations are used but the time series of mean daily CS and their 
corresponding mean XYZ values are obtained from station data to estimate a unique best. 
Thus, the slope of the MLR for the CS remained constant, but the y-intercept changes based 
on the mean CS and XYZ values. 


Trends in climate were calculated by comparing, through regression analysis, the annual 
mean CS in each county/subcounty against time (year) When the annual mean CS values 
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Table 1. 1960-99 Mean annual climate statistics and other independent variables used for this 
study, 
Climate statistic or other variable (and definition) Unit Variable name 
‘Mean annual precipitation intensity ea PRE 


{average for alt days) 


Mean annual precipration minus mean annual potential evapotranspration | sa meters per year | PRE-PET 
{average for all years} 


Mean annual minimum temperature i an 
{average for all days) Sonapi 

Mean annual mean temperature raises Grimis TME 
{average for all days) 

Mean annual maximum temperature dre clu TMX 
{average for all days) 

Mean annual number of wet days WDAY 
{days with measured precipitation) NS PEN 

Mean annual number of dry days ese ee DDAY 
{days without measured precipitation) ys pee rea 

Mean annual number of cold days Fae ac yee DAY 
(days with minimum temperatures below 0°C} ays pet ys 

Mean annual number of hot days 4 HDAY 
(days wih maximum temperatures above 35°C) al ions 

Mean annual sotar radiation oe R 
(average for all days) sg sili day, SAD; 

Mean annual total winter degree days ( {Tome = To} Where ‘ wo 
Tone = 3°C and Toe = {Tass + Ta) / 2. 2670 If negative de > 5 

Mean annual total summer degree days (4T me- Trae} Where “itn po 
Tran =3°C and Tue” {Tous + Toon) / 2, 2670 if negative) ANNE 

Mean elevation meters ELEV 

County area k AREA 
{total land area of county } ee (nes 


were missing or zero, no trend was calculated and a zero trend value was assigned to the 
county/subcounty Simulated CS in each county/subcounty for the years 1960 and 1999 were 
calculated using the trend regressions and the mean annual CS. Hence, the differences 
between the simulated CS values for the two years represent the magnitude of the trend over 
the 40 year time period and not the differences between any two years of actual CS data 


SOURCE AND PROCESSING OF ELEVATION AND AREA DATA 


Two additional variables used to augment the climate information for each 
county/subcounty were average elevation and total land area Elevation data were obtained 
from the USGS National Elevation Dataset [http fedc.usgs gov/products/elevation/ned html] 
and were projected from geographic coordinates referenced to the World Geodetic Survey of 
1984 to an Albers equal-area come projection using a bihnear mterpolation, 1000-meter cell 
resolution, and the following parameters ellipsoid — World Geodetic Survey of 1984, 1“ 
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standard parallel = 29 5°, 2° standard parallel — 45.5°, central meridian — -96.0°, latitude of 
origin = 23.0°, and no false easting or northing Average elevation was calculated as the mean 
of ali cells within each mappmg unit. Polygons for map units [http://www census. gov/geo/ 
wwwicob/scale html] were represented in an Albers equal-area conic projection using the 
same parameters as for the elevation data. Total mapping unit areas were determined from the 
county/subcounty polygons. 


STATISTICAL METHODS 


Multiple hnear regression (MLR) models (HELsEL & Hirscu, 1992) were developed 
using the SAS statistical software system (ANONYMOUS, 1990) to relate amphibian species 
richness to climate and location Dependent and independent model variables were standard- 
ized by subtracting the respective mean and dividing by the respective standard deviation. 
Standardized variables have equal weights in regression models, and the resulting model 
coefficients are proportional to their explanatory power in the models. “Best” and “stepwise” 
SAS regression procedures were used to screen potential models; however, neither method 
prevents correlated independent variables from entering the models. Multicollinearity among 
independent variables, as indicated by variance inflation factor (VIF) values greater than 10, 
can cause MLR model coefficients to be unrealistic in sign or magnitude (HELSEL & HIRSCH, 
1992) When an MLR model contained an independent variable with a VIF value greater than 
3, the independent variable was not used. 


Two sets of regression models were developed. one sct for the entire conterminous 
United States (including one model for anurans and one for urodeles), and one set for each of 
10 coarse-scale ecological regions (ANONYMOUS, 1997) (fig. 1). Because the primary objective 
of this research was to determine the degree to which climate explains patterns of amphibian 
species richness, model selection was manually supervised to favor climatic terms and prevent 
highly correlated independent variables from entering the same model. The adjusted coeffi- 
cient of determination (R°) and root mean square error (RMSE) statistics were used to 
evaluate the predictive skill of the models for a particular region or the nation (ANONYMOUS, 
1990; HetseL & Hirscu, 1992). The residuals between model and atlas estimates of species 
richness are used to compare the predictive capabilities of the national and regional models. 
Box plots are used to show the distributions of these residuals. The box plots show high and 
low outliers as circles. The central box extends from the 25" to 75" percentile of the data, and 
the box whiskers extend to the 5'" and 95™ percentiles. 


RESULTS 


NATIONAL REGRESSION MODELS 
Anuran species richness ranged from a maximum of 26 to a mimmum of 1 (fig 2a) The 


R’ for the national anuran model (fig 3a) was 0 64 (tab 2), with an RMSE of 3 07 species. 
Mean annual temperature and mean annual precipitation (fig 1b) accounted for the largest 
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Fig L Maps showing in the conterminous United States (a) coarse-level ecological regions and State, 
county and subcounty boundaries, and (b) 1960-1999 mean annual precipitation 
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Fig 2 Maps showing in the contermmous United States (a) anuran species richness and (b) urodele 
species richness, both from the Amphibian Research and Monstonng Intuative (AR MI} National 
Atlas for Amphibian Distributions, 
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Fig 3 - Maps showing in the conterminous United States the national regression model estimates of (a) 
dauran species richness and (b) urodele species richness. 
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Tab 2. Best-fitting national and ecological region standardized regression models of amphibian 
species nchness and adjusted coefficient of determination (2°) NA. North American: NW. 
Northwestern Sce tab. 1 for other abbreviations 


National or ecological Region Regression model R 
Ta | 
National 0 57*PRE + 0 56*TME - 0 49*PRE-PET + 0.07*ELEV 0.64 
Eastern Temperate Forests 0 71° FMX -0 21*ELEV +0.10*WDAY + 0 O7*AREA 063 
Great Plans 0 78*TME + 0 18*PRE + 0 08*AREA — 0 07*ELEV 078 
NA Deserts |0 79*TMX + 0 36*PRE + 0 22*AREA + 0.16*ELEV 047 
NW Forested Mountains -0 31* WDD - 0 27*ELEV + 0 25*AREA 023 
Norther Forests l-0 35*ELEV + 0 34*CDAY + 0 33*PRE 015 
Mediterranean California 0 45*ELEV - 0 41*CDAY + 0 39*RAD + 0 21*PRE 034 
Marine West Coast Forests 0 34*TME + 0.22*ELEV + 0 J4* AREA 025 
Temperate Sterras No statistically significant model. 3 
Southem Semi-Arid Highlands |0 80*ELEV + 0 68*RAD — 0 68*CDAY 068 
Tropical Wet Forests Too few mapping units (5) to develop a model. al 
Urodeles 

National 0.70*PRE - 0 24°PRE-PET + 0.11*TME - 0.05°ELEV 045 
Easter Temperate Forests  |-0.65"WDD + 032*ELEV + 0 29*WDAY + 0 15*PRE PLT 050 
Great Plans 0 55*PRE + 0.20* TME + 0 16*PRE-PLT + 0 15*AREA 049 
NA Deserts -0 57*ELEV + 0 37*PRE + 0 20*AREA - 0 19*TME 027 
NW Forested Mountains 055*PRE + 0 20*1ME -0 18*ELEV 0.60 
Northern Forests -0 67*WDD + 0 35*ELEV + 0 33*PRE on 
Mediterranean California O47*PRE-O43*WDD 0 34*TMX + 0 26*EI EV 050 
Marine West Coast Forests 0 33*PRE +0 31° TME + 0 23°LLEV + 0 22*RAD 040 
Temperate Sierras No statistically significant model. S 
Southern Semi-Arid Highlands | No statistically significant model. - 
Tropical Wet Forests Too few mapping umts (5) to develop a model a 


proportion of the variation (because they have the largest model coefficients, tab. 2), and were 
both positively associated with species richness Mean annual precipitation minus mean 
annual potential evapotranspiration also accounted for a substantial proportion of the 
variation and was inversely associated with species richness. Mean mapping unit elevation 
accounted for a small proportion of the variation and was positively associated with anuran 
species richness. The national regression modet overestimated anuran species richness along 
the Mississippi embayment and in parts of California, Florida. and Oregon The model 
underestimated anuran species richness along the Atlant. coastal plain and in parts of Mame 
and Texas (fig 2a, 3a). 


Urodele species richness ranged from a maximum of 30 to a minimum of 0 (fig. 2} The 
F for the national urodele model (hg 3b) was 0 45, with an RMSE of 4 54 species. The 
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national urodele model (tab 2) used the same CS as the national anuran model, but the 
coefficient values were appreciably different. Mean annual precipitation accounted for the 
largest proportion of the variation and was positively associated with species richness. Mean 
annual precipitation minus mean annual potential evapotranspiration accounted for a smaller 
proportion of the variation and was inversely associated with species richness, The national 
regressions model overestimated urodele species richness in the central United States and in 
parts of Florida and Washington; and underestimated species richness in most of the Eastern 
United States except for Florida and Maine (fig. 2b, 3b). 


REGIONAL REGRESSION MODELS 


Separate regression models were developed for each coarse-resolution ecological region 
(fig. 1a) to evaluate regional differences m the strength of chmate as a predictor of amphibian 
species richness. No models were developed for the Tropical Wet Forest ecological region, 
which was predominant only in five mapping units, and represented less than 0.3 % of the 
contermmous United States. The mean of anuran and urodele species richness for these five 
mapping units was used in place of a model. 


Eastern Temperate Forests 


The Eastern Temperate Forests ecological region was predommant in 1,789 mapping 
units, representing 31.8% of the contermsnous United States (fig. la) Anuran species 
richness ın this ecological region ranged from 2 to 26, and urodele species richness ranged 
from | 10 30 The R? for the Eastern Temperate Forest ecological region anuran model (fig. 4a) 
was 0.63, and the RMSE was 2.91 species. Mean annual maximum temperature accounted for 
the largest proportion of the variation and was positively associated with species richness, The 
Eastern Temperate Forest model overestimated anuran species richness in parts of Arkansas, 
Florida and Louisiana, and underestimated anuran species richness along the Atlantic 
Coastal Plain and in parts of Alabama and Indiana The residuals (model estimate minus 
atlas estimate} for the Eastern Temperate Forest anuran model were much smaller than the 
residuals between national model and atlas estimates of anuran species richness in those same 
mapping units (fig. 5a; gray box plots are residuals from regional models and black box plots 
are residuals from national model in the same mapping units) The R? for the Eastern 
Temperate Forest ecological region urodele model (fig 4b) was 0.50, and the RMSE was 3.78 
species. The total of mean annual winter degree days accounted for the largest proportion of 
the variation, and was inversely associated with species richness, The Eastern Temperate 
Forest model overesumated urodele species richness in parts of Arkansas, Florida, Illinois, 
Louisiana and Texas, and underestimated urodete species richness along the Eastern coastal 
and intand plains and in parts of Alabama, Indiana and Kentucky. The residuals for the 
Eastern Temperate Forest urodele model were smaller than the residuals between the national 
model and atlas estimates of urodele species richness (fig 5b). 
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Fig 4 Maps showing in the conterminous United States a compilation of reg.ona. regression models 
estimates of (a) anuran species richness and (b) urodele species richness. 
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GREAT PLAINS 


The Great Plains ecological region was predominant in 837 mapping units, representing 
28 9% of the conterminous United States (fig. La). Anuran species richness in this ecological 
region ranged from 1 to 23, and urodele species richness ranged from © to 13. The R’ for the 
Great Plains ecological region anuran model (fig. 4a) was 0 78, and the RMSE was 2.12 
species. Mean annual temperature accounted for the largest proportion of the variation and 
was positively associated with species richness. The Great Plains model overestimated anuran 
species richness in parts of Kansas, Missouri and Nebraska; and underestimated anuran 
species richness in parts of North Dakota, Oklahoma and Texas. The residuals for the Great 
Plains anuran model were much smaller than the residuals between the national model and 
atlas estimates of anuran species richness (fig Sa). The R° for the Great Plains ecological 
region urodele modet was 0.49, and the RMSE was 1.10 species. mean annual precipitation 
accounted for the largest proportion of the variation and was positively associated with 
species richness. The Great Plains model overestimated urodele species richness in parts of 
Towa, Kansas and Oklahoma, and underestimated urodele species richness in parts of North 
Dakota and Texas. The residuals for the Great Plains model were much smaller than the 
residuals between the national model and atlas estimates of urodele species richness (fig 5b) 


NORTH AMERICAN DESERTS 


The North American Deserts ecological region was predominant in 349 mapping units 
representing 19.8 ‘% of the conternunous United States (fig. la) Anuran species richness in 
this ecological region ranged from | to 16, and urodele species richness ranged from 0 to 4. 
The R’ for the North American Deserts ecological region anuran model (fig 4a) was 0,47, and 
the RMSE was 2 05 species. Mean annual maximum temperature accounted for the largest 
proportion of the variation and was positively associated with species richness. The North 
American Deserts model overestimated anuran species richness in parts of California and 
Utah, and underestimated anuran species richness in parts of Arizona and Texas, The 
residuals for the North American Deserts anuran model were much smaller than the residuals. 
between the national model and atlas estimates of anuran species nchness (fig 5a). The R° for 
the North American Deserts ecological region urodele model (fig 4b) was 0 27, and RMSE 
was 0 65 species Mean mapping unit elevation accounted for the largest proportion of the 
variation and was inversely associated with species richness. The North American Deserts 
model overestimated urodele species richness in parts of Cahfornia and Nevada, and under- 
estimated urodele species richness in parts of California. The residuals for the North Amen- 
can Deserts urodele model were smaller than the residuals between the national model and 
atlas estimates of urodele species richness (fig. Sb). 


NORTHWESTERN FORESTED MOUNTAINS 


The Northwestern Foresied Mountains ecological region was predominant in 289 
mapping units, representing 9 1 * «of the conterminous United States (fig la) Anuran species 
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richness in this ecological region ranged from | to 8, and urodele species richness ranged from 
Oto 12. The R? for the Northwestern Forested Mountains ecological region anuran model (fig. 
4a) was 0.23, and the RMSE was 1.20 species. The total of mean annual winter degree days 
accounted for the largest proportion of the variation and was inversely associated with species 
richness. The Northwestern Forested Mountains model overestimated anuran species rich- 
ness in parts of Colorado and Idaho, and underestimated anuran species richness in parts of 
Oregon. The residuals for the Northwestern Forested Mountains anuran model were much 
smaller than the residuals between the national model and atlas estimates of anuran species 
richness (fig. Sa) The R’ for the Northwestern Forested Mountains ecological region urodele 
model (fig. 4b) was 0.60, and the RMSE was 1.77 species. mean annual precipitation 
accounted for the largest proportion of the variation and was positively associated with 
species richness, The Northwestern Forested Mountains model overestimated urodele species 
richness in parts of Washington, and underestimated urodele species richness in parts of 
Colorado and Oregon. The residuals for the Northwestern Forested Mountains urodele 
model were much smaller than the residuals between the national model and atlas estimates of 
urodele species richness (fig. 5b). 


NORTHERN Forests 


The Northern Forests ecological region was predominant in 134 mapping units, repre- 
senting 5.2 % of the conterminous United States (fig. la) Anuran species richness in this 
ecological region ranged from 5 to 10, and urodele species richness ranged from 1 to 15. The 
R” for the Northern Forests ecological region anuran model was 0.15, and the RMSE was 0 92 
species. Mean mapping unit elevation accounted for the largest proportion of the variation 
and was inversely associated with species richness. The residuals for the Northern Forests 
anuran model were much smaller than the residuals between the national model and atlas 
estimates of anuran species richness (fig. 5a). The R° for the Northern Forests ecological 
region urodele model was 0.74, and the RMSE was 1.60 species. The total of mean annual 
winter degree days accounted for the largest proportion of the variation and was inversely 
associated with species richness. The residuals for the Northern Forests urodele model were 
stightly smaller than the residuals between the national model and atlas estimates of urodele 
species richness (fig. 5b). 


MEDITFRRANEAN CALIFORNIA 


The Mediterranean California ecological region was predommant in 277 mapping units, 
representing 2 1 "»of the conterminous United States (fig. La) Anuran species richness in this 
ecological region ranged from 2 to 9. and urodele species richness ranged from 0 to 10 The R? 
for the Mediterranean California ecological region anuran model was 0 34, and the RMSE 
was 1 04 species Mean mapping umt elevation accounted for the largest proportion of the 
vanation and was positively associated with species richness. The residuals for the Mediter- 
ranean California anuran model were much smaller than the residuals between the national 
model and atlas estimates of anuran species richness (fig. Sa} The R° for the Mediterranean 
California ecological region urodele model was 0 50. and the RMSE was 1.46 species. Mean 
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annual precipitation accounted for the largest proportion of the variation and was positively 
associated with species richness. The residuals for the Mediterranean Cahfornia urodele 
model were smaller than the residuals between the national model and atlas estimates of 
urodele species richness (fig. 5b). 


MARINE WEST COAST FORESTS 


The Marine West Coast Forests ecological region was predominant in 219 mapping 
units, representing 1.1 % of the conterminous United States (fig la) Anuran species richness 
in this ecological region ranged from 3 to 6, and urodele species nchness ranged from 2 to 10. 
The R for the Marine West Coast Forests ecological region anuran model was 0.25, and the 
RMSE was 0 67 species. Mean annual temperature accounted for the largest proportion of 
the variation and was positively associated with species richness. The residuals for the Marine 
West Coast Forests anuran model were much smaller than the residuals between the national 
model and atlas estimates of anuran species richness (fig. Sa). The R? for the Marine West 
Coast Forests ecological region urodcle model was 0 40, and the RMSE was 1 55 species. 
Mean annual precipitation accounted for the largest proportion of the variation and was 
positively associated with species richness. The residuals for the Marine West Coast Forests 
urodele model were much smaller than the residuals between the national model and atlas 
estimates of urodele species richness (fig. 5b). 


TEMPERATE SIERRAS 


The Temperate Sierras ecological region was predominant in 18 mapping units represen- 
ting 11°o of the conterminous United States (fig la). Anuran species richness in this 
ecological region ranged from 8 to 14, and urodele species richness was always 1. No 
statistically significant model of anuran species richness could be developed from the available 
independent variables (tab. 2). No model of urodele species richness was attempted since there 
was no variation in the dependent variable. 


SOUTHERN SEMI-ARID HIGHLANDS 


The Southern Senu-And Highlands ecological region was predominant in 17 mapping 
umts representing 0 6‘ » of the conterminous United States (fig 1a} Anuran species nchness 
in this ecological region ranged from 9 to 15, and urodele species richness was either 0 or 1 
The R for the Southern Semi-Arid Highlands ecological region anuran model was 0 68. and 
the RMSE was 0.88 species. Mean mapping unit elevation accounted for the largest propor- 
uon of the variation and was positively associated with species richness. The residuals for the 
Southern Senu-Arid Highlands anuran model were slightly larger than the residuals between 
the national model and atlas estimates of anuran species richness (fig 5a) No statistically 
significant model of urodele species richness could be developed from the available indepen- 
dent variables due to the limited variation in the dependent variable. 
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CLIMATE TRENDS 


Amphibian species richness was strongly associated with several of the mean annual 
climate variables, and mean annual precipitation and mean annual temperature were statisti- 
cally significant variables in 12 and 8 models, respectively (tab. 2) Increasing trends in annual 
mean temperature and precipitation were prevalent across much of the conterminous United 
States between 1960 and 1999 (fig. 6). Exceptions include decreasing mean annual tempera- 
ture in parts of the Great Plains ecological region, and decreasing mean annual precipitation 
in the southeastern part of the Eastern Temperate Forests ecological region. 


Discussion 


At the national level, the model for anurans performed better than that for urodeles (fig. 
5). Both models included mean annual precipitation as a strong variable for predicting 
patterns of richness. Kirstrr (1971) and DUELLMAN & Sweet (1999) previously noted a 
strong correlation in the conterminous United States between amphibian species richness and 
mean annual rainfall. Parutionmg the country by coarse-resolution ecological regions resul- 
ted in improved models for both anurans and urodeles. The residuals (model estimate minus 
atlas estimate) for the compilation of regional anuran and urodele models were much smaller 
than the residuals between national model and atlas estimates of anuran and urodele species 
richness in all mapping units (fig 5). In several cases the R? for the regional models were less 
than that of the national model, but the residuals were also smaller. In general, temperature 
variables (mean annual mean and mean annual maximum) figured more strongly in anuran 
models, whereas precipitation (mean annual precipitation intensity} had greater explanatory 
value in urodele models This makes sense from the perspective that there 1s no urodele 
counterpart to toads, hence, anurans are less restricted by amd conditions than are 
urodeles. 


In general, trends in chmate during 1960-1999 were toward wetter, warmer conditions for 
most of the conterminous United States. This could have provided more surface moisture 
availability for breeding habitat, and air and soil temperatures more amenable to regulating 
amphibian body temperatures throughout the year. Trends toward drier conditions in part of 
the southeastern United States and southwest Oregon may have resulted in reduced avatlabil- 
ity of breeding habitat in those areas. 


This effort to model the relations between anuran and urodele species richness and mean 
annual climate m the United States capitalized on the strong dependence of amphibians on 
their external environment for internal hydrothermal regulation. A limitation of the approach 
was that it assumed that the climate experienced by amphibians was reflected by long-term 
climate statistics summarized at the county/subcounty level. In fact, amphibians interact with 
climate at muluple scales, and alter their behaviors in concert with microhabitat features (sun 
flecks, burrows, duff. vegetation cover. wetlands. etc } to modify the effects of the broader- 
scale conditions. Therefore, the conditions represented by the data in this study likely 
addressed only the broadest effects of climate For this reason, the statistical models presented 
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here were aimed at the general level of anuran and urodele richness, and were not aimed at 
predicting the fate of particular species. 


Other explanatory variables may improve the ability to explain patterns of amphibian 
richness. For example, seasonal climate statistics may be more informative than annual 
statistics for certain measures, and additional landscape factors (e.g., current and historic and 
cover/use, hydrology, glaciation) and information such as evolutionary lineage could be very 
useful Additionally, models could be developed at the family level, or for groups of species 
having similar life history or developmental characteristics. 


The coarse-level ecological regions used for this study were somewhat problematic. 
Highly discontinuous mountain regions in the West often did not align well with 
county/subcounty units, so not all discontinuous portions of these regions were represented in 
the models. The largest regions (Eastern Temperate Forests, Great Plains and North Ameri- 
can Deserts) included a Jot of variation in temperature and moisture gradients. A finer level of 
regionalization (Level II regions defined in ANONYMOUS, 1997} may have been more appro- 
priate, as it would have subdivided the largest regions, while leaving the smaller regions intact, 
however, the number of map units per region may have been insufficient for developing 
models for several of the regions at this level. 


RÉSUMÉ 


Les amphibiens occupent une grande diversité d'habitats sur la planète, mais leur 
richesse spécifique est plus élevée dans les régions aux climats humides et chauds. Nous avons 
modélisé les relations statistiques entre la richesse spécifique en anoures et urodèles et le climat 
annuel des Etats Unis continentaux, et compare ces relations aux niveaux national et régional, 
Les variables modélisées ont été calculées pour des unités cartographiques correspondant aux 
contés ou aux sous-contés, et se sont appuyees sur des statistiques climatiques annuelles 
moyennes recueillies sur une période de 40 années (1960-1999), l'altitude moyenne et la 
surface des unités cartographiques, et des estimations de la richesse spécifique en anoures et 
urodeles, Les données climatiques ont éte obtenues à partir de plus de 7500 stations metéo- 
rologiques et ont eté mcorporées dans les données concernant les unités cartographiques au 
moyen de modèles de régression linéaire muluple Les richesses specifiques en anoures ct 
urodèles ont été calculées à partir de l'atlas national de distribution des amphibiens préparé 
par l'Amphibian Research and Monitoring Initiative (ARMI) de l'United States Geological 
Survey. Le modèle de régression hneare multivariée (MLR) national pour la richesse spéci- 
fique en anoures a un coefficient de déternunation ajuste (R°) de 0,64 et celui concernant les 
uradèles un R? de 0,45, Lorsque les Etats Unis sont divisés en regions ecologiques grossieres, 
on obtient des modèles pour les anoures dont les R` se répartissent entre 0,15 et 0.78 pour les 
anoures, et entre 0,27 et 0,74 pour les urodeles En général, les modèles régionaux pour les 
anoures se sont avéres plus fortement influencés par des variables de temperature, tandis 
que les variables liées à La precipitation avaient plus d'influence sur les modèles pour les 
urodeles, 
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